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In this work we investigate Ising and classical Heisenberg models for two and three dimensional 
lattices in presence of diluted ferromagnetic dimers. For such models the Curie temperature as a 
function of ratio of intra-dimer exchange coupling constant I a and other inter-site coupling constants 
Ib is calculated. In case dimer is treated exactly and environment within the mean-field approach it 
was found that even for Ia/Ib — > °o Tc remains finite. Similar analysis is proposed for rhombohedral 
phase of intermetallic compound Gd2Fei7 where so-called Fel-Fel "dumbbell" forms the diluted 
ferromagnetic dimer. It was shown that for such complex magnetic systems Tc is determined by all 
variety of exchange interactions and for the interval < Ia/Ib < oo Tc changes are not more than 
±10%. 

PACS numbers: 75.50.Ww, 78.20.-e, 71.20.-b 



I. INTRODUCTION 

There is a number of magnetic compounds with com- 
plex crystal structure and many magnetic ions per unit 
cell. Interatomic distances in such systems between near- 
est exchange-coupled ionic pairs could be quite different 
by its value. Thus the value of exchange coupling param- 
eter for a pair with shortest interatomic distance could 
be 2-3 times stronger than for a pair of atoms with larger 
interatomic distance. 

A typical example of such systems is R,2Fei7 class 
of intermetallic compounds. Here shortest Fel-Fel dis- 
tance corresponds to Fe ions in a so called "dumb- 
bell" positions (Wyckoff positions 4f for hexagonal 
ThzNii-j-type and 6c for rhombohedral ThzZnn-type 
crystal structures, see for details of magnetic struc- 
ture Ref. ll|). Based on band structure calculations 
within the LSDA+U method the exchange interac- 
tion parameter value for dumbbell pair in the hexag- 
onal phase Gd2Fei7 gives ferromagnetic exchange with 
7 U (1)=238.8 K for a distance r u (l)=2.400A^ While 
for Fe3-Fe3 pairs in the 12j Wyckoff positions for a 
distance of r33(l)=2.400 Acorresponding exchange pa- 
rameter value is found to be i33(l)=80.4 K£ Analo- 
gous calculations for rhombohedral phase of Gd2Fei7 
for Fel-Fel dumbbell pair gives ferromagnetic coupling 
with 7 n (l)=287.5 K (m(l)=2.385 A)£ For the next 
pair interatomic distance Fe2-Fe3 (9d and 18f Wychoff 
positions correspondingly) we obtained 723(1)=87.1 K 
(^23(1)=2.423 A)2. Thus one can see that coupling for 
dumbbell positions is 3 times stronger (in both phases) 
than an exchange coupling for the next smallest pair in- 
teratomic Fe-Fe distance. 

Let us note that these strong exchange bonds 7n(l) 
of Fel-Fel dumbbell positions are rather well spatially 
separated one from each other (rather diluted) and do 
not form any infinite magnetic cluster. Moreover these 
strongly coupled dumbbell pairs one can consider as some 
ferromagnetic dimers embedded into a infinite magnetic 



cluster formed by other weaker magnetic bonds. 

To this end there rises a question: how strong influ- 
ences presence of such diluted ferromagnetic dimers on a 
Curie temperature? Well could be that such dimers give 
dominating contribution to the Curie temperature of a 
complex magnetic system described above. 

Here we solve a problem of influence of diluted ferro- 
magnetic dimers on a Curie point. First for model lat- 
tices with Ising or classical Haisenberg spins with differ- 
ent numbers of lattice sites we vary exchange interaction 
strengths of the dimer from zero to infinity while all other 
exchange interactions remains finite (see Sec. |H]). Then 
we investigate rhombohedral phase of Gd2Fei7 using the 
same approach (Sec. IIII[) . Finally we conclude our work 
with a summary (Sec. IIV|) . 



II. THE CURIE TEMPERATURE FOR MODEL 
LATTICES WITH DILUTED FERROMAGNETIC 
DIMERS 

A. Ising spins case 

Let us consider square lattice with lattice parameter 
a. On each site we have Ising spin = ±1. We also 
specify that spins n and t 2 on a plaquette (selected on 
Fig. H] by dashed circle) form a dimer and are coupled by 
exchange integral I a >0 while each other pair of spins 
T i — T 3i r 3 — T 4 an d T4 — 72 are coupled with Ib >0 and 
Ib =h I A- Thus magnetic elementary cell contains four 
magnetic atoms and the lattice has period 2a. 

Using Weiss mean-field approximation it is easy to 
obtain four self-consistent equations for observable spin 
value o i —< Ti > (i=l,2,3,4) for each magnetic sublat- 
tice. 
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FIG. 1: Square lattice with spins in the lattice sites with the 
period 2a (magnetic unit cell surrounded by dashed circle). 
Strong exchange bond Ia — fat black line, thick lines — ex- 
change bonds Ib- Numbers 1,2,3,4 enumerate spins (magnetic 
sublattices) in the magnetic unit cell. 



cri = th{l3[I A a 2 +lB{(T2+2a 3 )}} 

o-2 = th{(3[I A ai+ I B {<Ji+2<j 4 )]} 

a 3 = th{/3I B (2a 1 +2a i )} 

g a = th{pi B (2a 2 +2a 3 )} (1) 



with /3 = l/k B T. 

If T —> Tc and <ti , a 2 , o"3 , 04 
linear equations determining Tc 



Eq. Q] reduces to 



k B Tc{lA = I B ) = 4I B . However in limiting case Ia —> 
00 Eq. (O leads to k B T c — > 00. 

It is obvious that the above mentioned result is not cor- 
rect from physical point of view. In case Ia S> I B strong 
ferromagnetic exchange coupling I a within a dimer pre- 
serves parallel alignment of spins n and t 2 . While de- 
coupling of weak exchange bonds I B (which form infi- 
nite ferromagnetic cluster) is responsible for destruction 
of ferromagnetic order. Thus even for I a — > 00 Curie 
temperature Tc should be finite and proportional to I B . 

The origin of such paradox is applied mean-field ap- 
proximation where strong fluctuations of exchange en- 
ergy because of spin flips are absent. The fluctuations be- 
come anomalously strong if we deal with groups of spins 
strongly coupled with each other. To this end one should 
use approach proposed by Scalapino, Imry and Pincus 3 
for quasi one dimensional magnetic systems with two sig- 
nificantly different exchange constants I a and I B . The 
main idea is to treat strongly coupled spins exactly while 
weak exchange interactions are treated within the mean- 
field approximation. 

For this purpose let us write down a Hamiltonian of 
such two spins cluster (dimer) in the elementary magnetic 
cell shown on Fig. [T] 



H c ; = -Iatit 2 - I b {<t 2 + 2ct 3 )ti 

- I B {(?1 + 2(Ta)T2. 



(4) 



Here interaction between spins t\ and r 2 in the dimer 
which are coupled with a large exchange constant I a is 
kept in the exact form and interaction of the same spins 
Ti and t 2 with others via small I B is treated within the 
mean-field approximation. 

Partition function for the dimer in the environment is 



Zd = 2{e 0lA ch[pI B {<n+ (T 2 4-2(T3 

+ e- f3lA ch[(3I B (a 1 -a 2 -2a 3 + 
Then for average values of spins a± 



f- 2(74)] - 
2(7 4 )]}. 



(5) 



=< Ti > and 

(72 =< t 2 > one can derive self-consistent set of equa- 
tions 



o-i = Pc[{Ia+I B )cti+2Ib(T3\ 

°i = f3 c [(lA + I B )cn+2I B a 4 } 

<rs = f3 c I B (2<r 1 +2<T 4 ) 

(74 = /3 C /s(2(72 + 2(73). (2) 

For ferromagnetic ordering in elementary magnetic cell 
shown on Fig. [1] follows relations o\ = a 2 and (73 = (74. 
Then 

k B T c = \{Ia + 3/b + ^l\-2I A I B + ni B ) (3) 

If I A = it comes out from Eq. © k B T c (I A = 0) = 
5(3 + Vl7)I B — 3.562/b and for Ia = I B one can find 



(71 = — le^S^/B^l +(7 2 + 2(7 3 + 2(74)] - 

- ((7! -(72 - 2(73 + 2(74)]}, (6) 

(72 = ^-{e l3lA sh[/3I B {<T 1 + (7 2 + 2(7 3 + 2(74)] + 
^cl \- 

+ e-P lA sh\pi B {ax- (7 2 - 2(7 3 + 2(7 4 )]}. (7) 

For average values of spins a 3 —< t 3 > and (74 =< T4 > 
weakly coupled to environment via I B still valid last two 
equations of Eqs. JT]). 

It is remarkable that for the limit I a — > 00 immediately 
follows 

(7l = (7 2 = th[pI B ((71 + (7 2 + 2(7 3 + 2(7 4 )] ■ (8) 



It means that for infinitely strong exchange coupling 
between spins 0\ —< t\ > and 02 —< r 2 > (in case 
the spins are strictly parallel to each other) environment 
mean- field I B (<J\ + 204) acting on spin 72 which belongs 
to the dimer is transferred onto spin t% and thus is added 
to Weiss field I B {u2 + 20-3). Similar "transfer" of Weiss 
field from spin t\ happens to spin t%. 

If we turn back to finite values of I a, considering 
(7i = 02 and 03 = 04 because of symmetry of magnetic 
surrounding, one can obtain two self-consisted equations 



0-3 



sh[l3I B {2a 1 +4g 3 )] 

-2j3I A 



ch[/3I B (2a 1 
th[2/3I B (tri - 



- 4ct 3 )] 
03)] • 



(9) 



In case of T —> Tc one can linearize Eqs. © for 
(Ti , 03 — > and get uniform linear set of equations 



<Tl = PclB[l + th{p C I A )]{<Tl 
0"3 = 2/3 C J B ((7 1 + 3 ). 



2a 3 ), 



(10) 



Since Eqs. (|10|) are valid only for I a ^> Ib, then for 
/a = 10/b one can find k B T c {I A = 10/ B ) = 4.792J B 
and for Ia — > 00 correspondingly k B Tc(jA — > 00) = 
2(1 + >/2)ii3 = 4.828I B . 

To summarize this subsection we have investigated a 
square lattice model with four spins in the elementary 
magnetic cell and one strong exchange bond Ia out of 
four. Varying the selected exchange coupling constant 
in the interval < Ia < 00 we obtained variation of 
the Curie point from k B Tc(lA = 0) = 3.562/b upto 
k B T c {I A -> 00) = 4.828/s. 



B. Elementary magnetic cell finite size effect on Tc 

Now one traces how an increase of number of mag- 
netic sites in the elementary magnetic cell changes Curie 
temperature in presense of the dimer. In case of square 
lattice with the periodicity of magnetic elementary cell 
equal to 3a containing nine magnetic sites (see Fig. [5]) 
magnetic structure is described by nine magnetic sublat- 
tices. However using symmetry of magnetic environment 
one can find following relations between different sublat- 
tice spin values o\ = 02, 04 = 05 = 07 = 08 and <j§ = 09. 
Thus if < Ia < Ib employing mean-field approxima- 
tion for T — > Tc one can obtain four linear equations 



01 = (3 c [Ia<Ji + is (03 + 204)], 

<?3 = fc[Is(2ai + 2a 6 )], 

04 = l3c[I B (cri+2a 4 + a 6 )}, 

°6 = /3 c [/b(03 + 20 4 + 6 )]. (11) 

By solving fourth order determinant one can get 
k B T c (I A =0) = 3.820/ B . 
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FIG. 2: The same as Fig. 1, but for the period 3a (9 magnetic 
sublattices) . 



In case of strong exchange bond I a ~> Ib it is necessary 
to write corresponding Hamiltonian similar to Eq. (U) 
using relations 04 = 05 = 07 = 0g 

H c /(1, 2) = Iati t 2 - ^5(03 + o- 4 + 0-7)7-1 - 

- Ib{<J3 + 05 + 08)72 = 

= -IaTit 2 - I B (a 3 + 20 4 )(ti + r 2 ). (12) 
Then average value 01 = < n > is 

sh[/3I B (2a 3 +Aa i )] 



01 



ch[pI B (2a 3 + Aa 4 )]+e-^ 



(13) 



In the vicinity of Tc Eq. (|13j) is linearized for small values 
0-3, 04 — ?> and gives 



ax = PcIb[1 + thtfcUWi + 20-4), 



(14) 



and should be used as a first equation in the Eqs. (|lll) . 
Then one can calculate k B Tc(lA — 10/_b) = 4.331/s and 
k B T c {I A -> 00) = 4.340/s. 

Similar calculations were done by us also for the case 
of two dimensional lattice with period 4a as well as for 
three dimensional models with elementary magnetic cells 
(2a x 2a x 2a) shown on Fig. [3] (8 magnetic sublattices), 
(3a x 3a x 2a) (18 magnetic sublattices) and (3a x 3a x 3a) 
(27 magnetic sublattices). Corresponding results are col- 
lected in Tabble HI Is is clearly seen that under increase 
of number of magnetic sites in the elementary magnetic 




Ising spins case 





Q— Q 2ax2a ( 4 atoms) 



□— □ 3ax3a ( 9 atoms) 
4ax4a (16 atoms) 

A-A 2ax2ax2a ( 8 atoms) 
3ax2ax2a (12 atoms) 

"^-^ 3ax3ax2a (18 atoms) 
3ax3ax3a (27 atoms) 
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FIG. 4: Dimensionless Curie temperature versus selected ex- 
change bond strength to other exchange bond strengths ratio 
Ia/Ib for two and three dimensional cubic lattices of Ising 
spins with different size of elementary magnetic cells. 



FIG. 3: The same as Fig. 1, but for the three dimensional 
cubic lattice (8 magnetic sublattices). 



cell and corresponding increase of number of weak mag- 
netic bonds influence of the dimer on Tc goes down and 
becomes insignificant even in the limit I a — > oo. 

Roughly speaking infinitely strong exchange bond 1a 
(in case spins are strictly parallel in the dimer) produces 
an effective spin with doubled value of magnetic moment. 
At the same time it decreases number of magnetic sites 
by one. After all such dimer contributes rather weak 
into free energy of the system in comparison with other 
spins. The same data is displayed on Fig. [4] Here one 
can see a tendency for the Tc to lesser and lesser deviate 
from Tc(Ia — Ib) for larger number of atoms in the 
elementary magnetic cell. 



C. Classical Heisenberg spins case 

In our days there is a way to compute exchange in- 
teraction parameters between different sites of Heisen- 
berg model with classical spins proposed by Lichtenstein 
et al. in Ref. HI- In this approach exchange parameters 
are calculated as a second derivative of total energy with 
respect to small angles of magnetic moments deviation 
from collinear magnetic configuration. Thus one should 
understand how Curie temperature of ferromagnetic lat- 
tice of classical Heisenberg spins depends on presence of 
diluted magnetic dimers. Here we present analysis simi- 
lar to Sec. Ill Al 

At the beginning we consider a lattice shown 
on Fig. [TJ with classical spin vectors Si = 
S(sindiCOS(pi,sin9iSin(j)i,cos9i). Now we calculate 
average values of spin vectors for ferromagnetic case 
assuming that spontaneous magnetic moment is directed 



along z-axis. It means that < Sf >—< Sf >— and 
mi =< Sf >= S < cosOi >^ 0. Within the mean-field 
approximation for spins average value on sites of each 
of four magnetic sublattices one can get following 
equations: 



mi = SC(f3S[lAm 2 + I B {m 2 + 2m 3 )]), 

m 2 — S£((3S[lAmi + lB(mi +2m4)]), 

m 3 = SC(pSI B [2mi + 2m i ]), 

m 4 = SC(/3SI B [2m 2 + 2m 3 ]). (15) 



where C(x) — cth{x) 
T c C{x) « f 



— - Langevin function. For T 



mi 



m 2 



m 3 



!7i4 = 



x <C 1 (compare with Eq. [5]): 
- Pc y [(Ia 



lB)m 2 + 2/37713], 



p c — \[I A + I B )mi + 2I B m 4 ], 
S 2 

/3cy/s(2m 1 + 2m 4 ), 
S 2 

/3 c —I B (2m 2 + 2m 3 ). 



(16) 



Exploiting symmetry relations mi = m 2 and 777,3 = TO4 
one can end up with equation on Tq 



k B T c = ~(h 



3J S 



2I A I B + 17I 2 B ), (17) 



n2 

which differs from Eq. ([3| only with coefficient ^- on the 
left side. Similar to Ising case simply mean-field treat- 
ment gives incorrect results, for the case — > 00 namely, 
Tc — > 00. 

Let us employ once more the approach of Ref. 3, where 
for the case I a S> I B for a dimer we write down a Hamil- 
tonian with exchange interaction in the dimer is taken 
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TABLE I: Curie temperature value for different magnetic elementary cells with Ising spins with selected strong exchange bond 
I a as a function of ratio Ia/Ib- Tc is given in units Is- 



Magnetic 


Magnetic 












elementary 


atoms 


Za=0 


Ia=Ib 


Ia=5I b 


Ia=WIb 


I A — >■ OO 


cen 


number 












la x 2a 


4 


3.562 


4 


4.582 


4.792 


4.828 


3a x 3a 


9 


3.820 


4 


4.254 


4.331 


4.340 


4a x 4a 


16 


3.905 


4 


4.149 


4.193 


4.199 


2a x 2a x 2a 


8 


5.785 


6 


6.496 


6.785 


6.899 


3a x 2a x 2a 


12 


5.862 


6 


6.350 


6.562 


6.641 


3a x 3a x 2a 


18 


5.910 


6 


6.241 


6.389 


6.441 


3a x 3a x 3a 


27 


5.929 


6 


6.196 


6.319 


6.362 



exactly. The exchange coupling constant is I a and spin 
vectors S\ S 2 ■ Coupling of these two spins to other ones 
via weak exchange integral Ig we consider within mean- 
field approximation: 

H d (l, 2) = -I A Sx§ 2 - h x {Sl + SI), (18) 

here Weiss field is hi = /g(mi + 2ms). Corresponding 
partition function is 



Z cl (l, 2) = / dfii Jdn 2 exp{^[I A S 1 S 2 + hi(Sf + SI)}} = 

— 4n J dQiexp((3Shicos8i) x 

sh{pS^I\S 2 + 2I A Sh 1 cos9 1 ~Th 1 1 ) 
X PS^JP A S 2 + 2I A Shicos8i + h'{ ' 

where dQi — sin9id9id(j)i is solid angle element with < 

6i < 7T, < (j>i < 27T. 

Since to get Tq one needs to linearize equations for m\ 
and 7713 (that is hi). Therefore we calculate Z c i(l,2) up 
to the order h\. 

+ l -[i + c(pi A s 2 wshiY + ... }. (20) 

Then average value of spin mi on site 1 in the elementary 
magnetic cell is 

mi = i(mi + m 2 ) = i(< 5 X Z > + < >) = 



2 

1 1 



9 



2^(1,2)9(^1 



-^(1,2). 



After differentiating of Eq. (f20| and following expanding 
of self-consistent equation for 1713 around small values of 
mi and 1TI3 we come to two equations 



mi - /3c^-3-[l + A/3c^5 2 )](mi+2m 3 ), 
m 3 2/5 c ^B-^-(mi + m 3 ). 



1.5- 



Heisenberg spins case 



O-O 2ax2a ( 4 atoms) 
Q— □ 3ax3a ( 9 atoms) 
<5 — ^> 4ax4a (16 atoms) 
A-A 2ax2ax2a ( 8 atoms) 
3ax2ax2a (12 atoms) 
3ax3ax2a ( 1 8 atoms) 



3ax3ax3a (27 atoms) 




Infinity 



FIG. 5: The same as Fig. 4, but for classical Heisenberg spins. 



(21) 



Last equations for Tc coincides with Eqs. (|10|) for Ising 

spins upto substitutions Ib^ > Ib and £((3cIaS 2 ) — > 

th(PclA)- It gives us a way to transform correspond- 
ing equations determining Tc from Sec. Ill Al to the 
case of classical Hesenberg spins in frame of the same 
model approach. Let us also mention that limits 
lim/^oo C(/3 c IaS 2 ) = 1 and linij^oo th(/3 c lA) = 1 co- 
incide in magnitude, while for finite values of I a values of 
MPcIa) and C{(3 c IaS 2 ) are slightly different (if S = 1 
is assumed). In the Table ILT1 values of Tc are listed for the 
case of S = 1 for different strength of exchange coupling 
within dimer I A - From Fig. [5] one can see that values 
of Tc vary not more than 10% from Tc(Ia = Ib) when 
< I a < oo. Again with increase of number of atoms in 
elementary magnetic cell Tc is getting closer and closer 
to T c (I a =Ib). 



III. INFLUENCE OF DUMBBELL EXCHANGE 
STRENGTH ON T c OF THE RHOMBOHEDRAL 
PHASE OF GD2FE17 

Here we examine influence of strongly exchange cou- 
pled dumbbell Fe atoms on Tc of the rhombohedral phase 
of Gd2Fei7. In this compound there are 17 magnetic 
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TABLE II: Curie temperature value for different magnetic elementary cells with classical Heisenberg spins with selected strong 
exchange bond I a as a function of ratio Ia/Ib for S = 1. Tc is given in units Ib- 



Magnetic 


Magnetic 












elementary 


atoms 


Ia=0 


Ia=Ib 


Ia=5Ib 


Ia=WIb 


I A — >■ CO 


cell 


11U111UC1 












la x 2a 


4 


1.187 


1.333 


1.489 


1.547 


1.609 


3a x 3a 


9 


1.273 


1.333 


1.389 


1.407 


1.447 


4a x 4a 


16 


1.302 


1.333 


1.372 


1.386 


1.400 


2a x 2a x 2a 


8 


1.928 


2 


2.144 


2.215 


2.299 


3a x 2a x 2a 


12 


1.954 


2 


2.101 


2.152 


2.214 


3a x 3a x 2a 


18 


1.970 


2 


2.068 


2.104 


2.147 


3a x 3a x 3a 


27 


1.976 


2 


2.056 


2.085 


2.121 



(Fe) atoms in the elementary magnetic cell, which oc- 
cupies four types of nonequivalent crystallographic po- 
sitions. Correspondingly they have different surround- 
ings of neighboring Fe atoms. Namely, for Th2Znn~ 
type rhombohedral structure there are 2 Fel atoms in 
6c Wychoff positions (dumbbell positions), 3 Fe2 atoms 
(9d), 6 Fe3 atoms (18f) and 6 Fe4 atoms (18h). Lo- 
cal magnetic moments of different classes of Fe atoms 
are a bit different from each other: fiFei = 2.19/U.b, 

llFel = 2.26(13, MFe3 = 2.Y1[L B and /i Fe4 = 2.31/iB, 2 

Exchange interaction parameters for the first coordina- 
tion sphere of different Fe atom classes calculated within 
the LSDA+U approach^ in Ref. [H are presented in Ta- 
ble Iffll 



It was found that largest parameter of exchange in- 
teraction iii(l) = 287.5 K couples two Fel atoms 
in the dumbbell position with a distance inbetween 
rn=2.385 A. Next largest exchange coupling parameter 
J 4 4(l) = 182.2 K couples magnetic moments of Fe4 atoms 
with a distance r44(l)=2.490 A. Having at hand values 
of exchange interaction constants and local magnetic mo- 
ments for all Fe magnetic sublattices it is straightforward 
to estimate Tc using nearest neighbors mean-field ap- 
proximation. For this purposes absolute value of spin 
vector is Si = fip e ^/2fiB (i=l,2,3,4). To do that first 
one should define a set of self-consistent equations for 
average value of z-component of the spin m; =< Sf >, 
which are 



\k B TJ 



(22) 



where hi is Weiss field acting on spin Si from nearest 
neighbors. Next one should linearize right hand sides of 
Eqs. [22] expanding Langevin function for T — > Tc 



Sl_hi_ 
3 k B T' 



(23) 



Taking into account rhombohedral crystal structure of 
Gd2Fei7 one can derive set of linear equations for rrii 



defining Tc 

S 2 

mi = /3c— 
S 2 

m 2 =fcy 



S 2 

m 3 = /3c 



m 4 = /3c 



Sj 



[ In(l)«n(l)mi + Iia(l)«ia(l)ma 

+ Ii3(l)«is(l)m 3 + Iii(l)zu{l)mi} 

[ /21 (l)«2i (l)mi + 123(1)^23(1)^2 

+ ^33(1)233(1)^3], 

[ 131(1)231(1)^-1 + ^32(1)232(1)^2 

+ 133(1)233(1)^3 + 

+ (^34(1)234(1) + / 3 4(2)z 3 4(2))m 4 ], 

[ hi (1)^41 (l)wi + 142(1)242(1)^-2 

+ {h 3 (1)«43 (l)+^43 (1)^43 (l))m 3 4 

+ /44(1)244(1)™4], 



(24) 



corresponding exchange integrals hjiX) between Fei and 
Fej crystallographic classes in the first coordination 
sphere and nearest neighbors number 2^(1) on a distance 
rij(l) are taken from Table [TTT1 At that for cumbersome 
crystal structures 2^(1) ^ Zji(l) but ni^y(l) = rijZji(l), 
where rij number of atoms of sort Fei in a elementary 
magnetic cell. 

As was obtained in Ref. H Eqs. ((Ml) give T c =429 K 
for rhombohedral Gd2Fei7 which is slightly smaller than 
experimental one 475 Ki. A reason of this discrepancy is 
mainly absence in our model of exchange interactions be- 
tween Fe and Gd sublattices and oscillating exchange in- 
teractions with next coordination spheres. However this 
result shows that largest part (leading contribution) of 
Tc comes from exchange interactions of Fe sublattice be- 
tween nearest neighbors. Thus it is interesting to explore 
influence of strongest exchange /n(l) "in dumbbell" on 

At the beginning hypothctically one switches off the 
interaction Jii(l). Then solution of Eqs. ([2"4"|) leads to 
T c (/ii(l) = 0)=414 K which is only 3.5% less than in 
Ref. |2j. Now we consider opposite case hi(l) ^S> ^ij(l), 
i.e. consider dumbbell pair as dimer. Thus one can apply 
approach described in Sec. [Til Fel-Fel cluster with strong 
exchange coupling -Tii(l) will be treated exactly while 
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TABLE III: Parameters of exchange in the rhombohedral structure of GcbFeiT for the ions of the first coordination sphere. 



N 


Exchange (K) 


Distance (A) 


Number of neighbors 


Type 


1 


7n(l)=287.5 


* ii i . 




2ll(l) = l 


Fel (dumbbell) - Fel (dumbbell) 


2 


7 44 (1)=182.2 




=2.490 


Z 44 (l)=2 


Fe4 (corrugated plane) - Fe4 (corrugated plane) 


3 


J 34 (l) = 125.9 


r 34 (l) 


=2.551 


« 34 (1)=«4 3 (1)=2 


Fe3 - Fe4 (corrugated plane) 


4 


7 24 (1)=121.0 


r 24 (l) 


=2.448 


224(1)=4, 2 42 (1)=2 


Fe2 (corrugated layer) - Fe4 (corrugated layer) 


5 


I 3 4(2)=105.7 


r 34 (2) 


=2.613 


z 34 (2)=^43(2)=2 


Fe3 - Fe4 (corrugated layer) 


6 


7i 4 (l)=88.2 


r 14 (l) 


=2.639 


Zl 4 (l)=3, 2 4 l(l) = l 


Fel (dumbbell) - Fe4 (corrugated layer) 


7 


7 23 (1)=87.1 


r-23(l) 


=2.423 


223(1)=4, Z 32 (l)=2 


Fe3 (dumbbell) - Fe2 (corrugated layer) 


8 


7 i2 (l)=83.6 


ri 2 (l) 


=2.602 


2l 2 (l)=3, 2 2 l(l)=2 


Fel (dumbbell) - Fe2 (corrugated layer) 


9 


7i 3 (l)=74.1 


ri3(l) 


=2.740 


2 13 (1)=6, 2 3 l(l)=2 


Fel (dumbbell) - Fe3 


10 


7 33 (l)=-36.5 


r-33(l) 


=2.466 


z 33 (l)=2 


Fe3 - Fe3 



other exchange bonds within mean-field approach. After 
usual linearization of equations on m\ =< S{ > for small 
values of mi (i = 1,2,3,4) close to Tq one gets equation 

mi = Pc^[l+C{Pchi{l)Sl)\ x (25) 

x [7"i2 (l>i2 (l)m 2 + 7"i 3 (l)^i3(l)m 3 + 
+ 7i 4 (l)zi 4 (l)m 4 ], 

which should substitute the first equation of Eqs. (|24|) . 

For the rhombohedral phase of Gd2Fei7 second 
largest exchange parameter is I 44 (l)=182.2 K. Let us 
take 7n(l) = 10/ 44 (1)=182 2 K then T c {hi{l) = 
10J 44 (1))=456 K. And finally for the limit 7 U (1) ->■ 
oo, when C(/3cIn(l)S 2 ) —> 1, Curie temperature is 
Ib(iii(l) -> oo)=464 K. Thereby even dumbbell ex- 
change in(l) — > oo gives rise of Tc only 10%. Also from 
Eqs. (|26l) one follows that for iii(l) — > oo between spins 
Si and S2 Weiss field acting from other spins of atoms 
Fc2, Fe3 and Fe4 is doubled. 

Another interesting observation one can do for Curie 
point Tc calculated from averaged over all types of Fe 
atoms values of exchange constant I, spin value S and 
number of nearest neighbors z 

k B T c = ll(S) 2 z, (26) 

4 4 

S = ^n^/^n, = 1.12, 

i=i i=i 

4 4 

Z = ^ n i Z ijl X!" 1 = 10 ' 
i=l,j i=l 

I = ik B T c /{S) 2 z = 102.6 K, 

here K. From Table [TTTI one can see that most of 

values of exchange integrals are very close to the value 
of I. Thus these Fe-Fe exchange bonds are responsible 
for the Curie point value rather than Fel-Fel exchange 
bond only, independently how strong it is. 



IV. CONCLUSION 

In this work we proposed description of several model 
magnetic structures with diluted magnetic dimers which 
do not form any infinite magnetic cluster. The selected 
strong exchange bond I a is supposed to be much larger 
than other exchange interactions Ib in the elementary 
magnetic cell. Following ideas of Ref. [H the magnetic 
dimer is treated exactly while other couplings are treated 
within mean-field approximation. In contrast to regular 
mean-field approximation latter approach allows one to 
obtain finite Curie temperatures Tc even for I a — > 00 
which is physically correct. Also if was demonstrated 
that for < Ia/Ib < 00 Tc deviates just about ±10% 
from the value Tc(Ia = Ib) and is getting closer to that 
while number of atoms in the elementary magnetic cell 
grows. For the case 1a 00 we obtained doubling of 
spin value of spins forming a dimer, which corresponds 
to doubling of Weiss field acting on the spin from other 
spins in the elementary magnetic cell. After all one can 
conclude that such diluted magnetic dimers do not influ- 
ence very much on Tc value of the whole system in the 
case I a — > 00, and then even less for finite Ia- 

Based on these results we perform analogous calcula- 
tions of Tc for real system - the rhombohedral phase 
of Gd2Fei7. There is so called dumbbell Fel-Fel dimer 
with the largest exchange interaction value 7ii(l) in the 
system. We showed in this investigation that for such 
complicated crystal structure Tc is mainly defined by 
weaker exchange interactions of Fel with other Fe atoms 
in the elementary magnetic cell rather than by the 7n(l) 
exchange only. 
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